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Abstract. We report results from 2 1/2 yr of high precision radial velocity (RV) monitoring of Barnard's star. 
The high RV measurement precision of the VLT-UT2+UVES of 2.65 m s _1 made the following findings possible. 

(1) The first detection of the change in the RV of a star caused by its space motion (RV secular acceleration). 

(2) An anti-correlation of the measured RV with the strength of the filling-in of the H a line by emission. (3) Very 
stringent mass upper limits to planetary companions. Using only data from the first 2 years, we obtain a best-fit 
value for the RV secular acceleration of 5.15±0.89 m s _1 yr -1 . This agrees within 0.95<r with the predicted value 
of 4.50 m s _1 yr -1 based on the Hipparcos proper motion and parallax combined with the known absolute radial 
velocity of the star. When the RV data of the last half-year are added the best-fit slope is strongly reduced to 
2.97 ± 0.51 m s _1 yr -1 (3.0a away from the predicted value), clearly suggesting the presence of additional RV 
variability in the star. Part of it can be attributed to stellar activity as we demonstrate by correlating the residual 
RVs with an index that describes the filling-in of the H a line by emission. A correlation coefficient of —0.50 indicates 
that the appearance of active regions causes a blueshift of photospheric absorption lines. Assuming that active 
regions basically inhibit convection we discuss the possibility that the fundamental (inactive) convection pattern in 
this M4V star produces a convective redshift which would indicate that the majority of the absorption lines relevant 
for our RV measurements is formed in a region of convective overshoot. This interpretation could possibly extend 
a trend indicated in the behaviour of earlier spectral types that exhibit convective blueshift, but with decreasing 
line asymmetries and blueshifts as one goes from G to K dwarfs. Based on this assumption, we estimate that 
the variation of the visible plage coverage is about 20%. We also determine upper limits to the projected mass 
m sin i and to the true mass m of hypothetical planetary companions in circular orbits. For the separation range 
0.017 — 0.98 AU we exclude any planet with msini > 0.12 Mj up iter and m > 0.86 Mj up iter. Throughout the 
habitable zone around Barnard's star, i.e. 0.034 — 0.082 AU, we exclude planets with msini > 7.5 MEarth and 

TTl > 3.1 MNcptuno- 
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* Based on observations collected at the European Southern 
Observatory, Paranal, Chile (ESO programmes 65.L-0428, 
66.C-0446, 267.C-5700, 68.C-0415, and 69.C-0722) 



1. Introduction 

The search for extrasolar planets via precise stellar radial 
velocity (RV) measurements has led to more than 100 dis- 
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covery announcements up to now (e.g. Mayor & Queloz 
1995; Marcy & Butler 1996; Kiirster et al. 2000). A good 
fraction of these discoveries is summarized in the paper by 
Butler et al. (2002). The majority of the extrasolar planets 
has so far been found around main sequence stars of spec- 
tral types from late-F through K. There is only one M 
dwarf with planetary companions known, GJ876, which 
is orbited by two Jovian-type planets with orbital peri- 
ods of 30.1 and 61.0 day, respectively (Marcy et al. 2001; 
Dclfosse et al. 1998). Combining RV data with new astro- 
metric measurements Benedict et al. (2002) determine a 
mass of 1.89 ± 0.34 Mj up i te r for the outer planet. 

It is due to the intrinsic faintness of M dwarfs that 
this most abundant type of star tends to be underrepre- 
scnted in RV searches as very high measurement precision 
requires large telescopes. This introduces a strong bias in 
any attempt to estimate the frequency of planetary com- 
panions to stars as a whole. 

On the other hand, due to the low mass of M dwarfs 
an orbiting planet of a given mass causes a higher reflex 
motion of the star. Therefore, the present day measure- 
ment precision for differential radial velocities (DRV) 1 of 
a few m s _1 holds the prospect to discover terrestrial plan- 
ets (with a few Earth masses) in close-in orbits around 
the lower mass M dwarfs, e.g. in their habitable zones 
(Schneider 1998; Endl et al. 2003). 

For instance, for a main sequence star with a mass 
of 0.16 M Q , such as that of Barnard's star, the subject 
of the present paper (mass from Henry et al. 1999), the 
habitable zone lies at 0.034 — 0.082 AU (see Kasting et 
al. 1991) corresponding to orbital periods in the range 
5.75 — 21.5 d. Around a star of this mass a planet in a 
circular orbit with a projected mass msinz = 1 MEarth 
would cause a peak-to-peak RV amplitude of 2.4 m s _1 at 
the inner edge of the habitable zone and 1.6 m s' 1 at its 
outer edge. For msini = 3 Msarth the resulting peak-to- 
peak RV amplitudes would be 7.3 m s _1 and 4.7 m s^ 1 , 
respectively. 

When attempting to measure such small RV signals it 
becomes important to consider other effects that affect the 
determination of DRVs. Two such effects are important 
for M dwarfs. First, the high proper motion of nearby 
stars causes a secular change in their RV which for some 
stars can reach several m s _1 within a few years. Second, 
intrinsic variability due to stellar activity can mimic RV 
variations. Little is known so far about the level at which 
this complicates DRV measurements in M dwarfs. 

This paper reports on high precision DRV measure- 
ments obtained with VLT-UT2+UVES of Barnard's star 
which, as we demonstrate, exhibits a considerable secular 
change of its RV. It is known as a quite inactive star and 
was selected with the expectation that its DRV measure- 
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1 Throughout this paper we will differentiate between ab- 
solute radial velocity (ARV) and differential radial velocity 
(DRV), the latter being relative to any one of the measure- 
ments. It is the DRV for which high precision can be achieved 
and only DRVs are needed when RV changes are to be studied. 
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Fig. 1. Lower panel: Geometry of the space motion of a pass- 
ing star. v t and v T are, respectively, the transversal and radial 
components of the space velocity v. The current distance is d, 
the minimum distance is d Q , and vt is the trajectory from the 
current position to the point of nearest approach. Upper panel: 
The ARV as a function of time and its asymptotic limits ±w. 

ments are largely unaffected by activity; as we shall see 
this is not the case. 

The paper is structured as follows. We summarize the 
known properties of Barnard's star (Sect. 2) and proceed 
with a description of the concept of RV secular accelera- 
tion (Sect. 3). Our observations are described in Sect. 4 
along with a brief outline of our data modelling technique. 
In Sect. 5 we present results on the RV secular accelera- 
tion, on the relation between DRVs and stellar activity as 
evidenced in the strength of the H Q line, and on period 
analysis in the RV and H a data; we also provide mass up- 
per limits to planetary companions. We discuss our results 
in Sect. 6 and summarize our conclusions in Sect. 7. 

2. Properties of Barnard's star 

The M4V dwarf Barnard's star (GJ 699, HIP 87937) is the 
nearest star after Proxima Centauri and the a Centauri 
binary. Its parallax has been determined by Hipparcos 
as 549.01 ± 1.58 mas (ESA 1997; Perryman et al. 1997), 
i.e. the distance is 1.8215 ± 0.0052 pc. Barnard's star 
has long been known as the highest proper motion star 
in the sky (Barnard 1916, 1917). For epoch 1991.25 
Hipparcos gives \i a = —797.84 mas yr _1 and [is = 
10326.93 mas yr _1 for the proper motion in right ascen- 
sion and declination, respectively (ESA 1997). The tan- 
gential component of the space velocity (perpendicular to 
the line of sight) is thus Vt — 89.440 km s _1 . The secular 
change of the proper motion is 1.2 mas yr~ 2 (Gatewood 
1995). 

Barnard's star also has a high ARV. Dravins et al. 
(1999) find an astrometrically derived ARV of v r = 
— 101.9±6.5 km s~\ somewhat smaller (in absolute value) 
than the spectroscopic value of —110.506 km s _1 by 
Nidever et al. (2001) who estimate their systematic errors 
in their treatment of gravitational redshift and convective 
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Table 1. Kinematic data for Barnard's star. For all data used 
in our calculation (Hip = Hipparcos, Nid = Nidever et al. 2001, 
UV = our UVES data) we list the epoch, the distance d, the 
ARV v r , and the secular acceleration of the RV dv T /dt. We also 
compare our results for the point of closest approach (CA) with 
those of Garcia-Sanchez et al. 2001 (GS). 



Data/ 
Event 


Epoch 


d 

[PC] 


Vr 

[kins" 1 ] 


dv r /dt 
[ms^yr -1 ] 


Hip 


1991.25 


1.8215 


-110.532 


4.491 


Nid 


1996.91 


1.8209 


-110.506 


4.496 


UV 


2001.51 


1.8203 


-110.485 


4.499 


CA 
GS 


11730 
11700±ioo 


1.1458 

1.144±0.005 


0.0 
(0.0) 


18.043 
not given 



blueshift) to 0.4 km s . Other spectroscopic literature 
values are, e.g., — 110.9±0.2 km s _1 (Garcia-Sanchez et al. 
2001) and -109.7 km s _1 (Turon et al. 1998). 2 Thus the 
overall space velocity is v s=s 142 km s _1 . In the following 
we will adopt the ARV value by Nidever et al. (2001) not- 
ing that none of the results presented in this paper would 
be significantly altered, if, e.g., we took the astrometric 
ARV value by Dravins et al. (1999). The same would be 
true, if we replaced the Hipparcos astrometric data by the 
recently refined values by Benedict et al. (1999) who used 
the Fine Guidance Sensor 3 of the Hubble Space Telescope 
for interferometric fringe tracking astrometry. 

Due to its high space motion and sub-solar metallicity 
Gizis (1997) classifies Barnard's star as an "Intermediate 
Pop. II star" (following Majewski 1993), i.e. an object be- 
tween a halo star and a thin disk star. Its low observed 
X-ray luminosity (Marino et al. 2000; Hiinsch et al. 1999) 
and its presumable rotation period of 130.4 d (from weak 
evidence for w 0^01 photometric variations by Benedict 
et al. 1998) also indicate an inactive old star. 

For many years Barnard's star was also a prime candi- 
date for displaying astrometric pertubations attributable 
to one or two planetary companions in long-period or- 
bits of about 12 and 20 yr (van dc Kamp 1963, 1982). 
This was later refuted by Gatewood (1995). The most re- 
cent astrometric measurements by Benedict et al. (1999) 
placed stringent mass upper limits of 2.1 Mj up it r down to 
0.37 Mj U pit er for orbital periods of 50 d up to 600 d (see 
Fig. 8 in Sect. 5.4). 

3. The secular acceleration of the radial velocity 

The proximity and high space velocity of Barnard's star 
implies a detectable secular change of its radial velocity 
(RV secular acceleration; van de Kamp 1977). To illustrate 
this effect Fig. 1 (lower panel) shows the trajectory past 
the barycentcr of the solar system of a star in linear space 
motion, i.e. ignoring the subtle variations of the galactic 



2 The Hipparcos and Tycho Catalogues give the ARV value 
used to correct for the astrometric perspective acceleration 
(change in proper motion) of Barnard's star as —111.0 km s _1 . 



potential (see Garcia-Sanchez et al. 2001). The minimum 
distance d at the point of closest approach is given by 
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(i) 



where d, v t , and w r are the distance, and the tangential and 
radial velocity components, respectively, at some earlier 
epoch. The time t at the earlier epoch is 



t = t - 



where v = ^/v 2 4 
time at closest approach. Setting t, 
function of time is given by 
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is the space velocity and t Q is the 
= the ARV as a 
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The upper panel of Fig. 1 shows a graph of this function. 
Differentiation yields the secular acceleration of the RV: 



dv T (t) 
dt 



vH 2 



^JvH 2 + dl (vH 2 + <P )V* ■ 



(4) 



Kinematic data for Barnard's star using Eqs. (1-4) with 
the input data from Sect. 3 are summarized in Table 1 that 
lists the distance, the ARV, and its secular acceleration 
for the individual epochs of all measurements relevant to 
this paper. 3 For the time of closest approach to the solar 
system our results are compared with those of Garcia- 
Sanchez (2001). 

At the mean time of our UVES observations pre- 
sented here (epoch 2001.51) the RV secular acceleration 
is 4.499 m s _1 yr -1 . Over the time baseline for our UVES 
data of 2.422 yr the acceleration value increases by only 
0.0020 m s -1 yr -1 ; we neglect this small increase and 
treat the secular acceleration as a purely linear function 
when comparing it to our UVES data. 

4. Observations and data modelling 

Barnard's star was observed with the VLT-UT2+UVES 
in the framework of our ongoing precision RV survey of 
M dwarfs in search for extrasolar planets. UVES was used 
together with its self-calibrating iodine gas absorption cell 
operated at a temperature of 70° C. Image slicer #3 and 
an 0.3" slit were chosen yielding a resolving power of R = 
100 000 - 120 000. At the selected central wavelength of 
600 nm the useful spectral range containing iodine (I 2 ) 
absorption lines (~ 500 — 600 nm) falls entirely on the 
better quality CCD of the mosaic of two 4 K x 2 K CCDs. 

Our analysis of the stellar spectra to determine differ- 
ential radial velocities (DRV) followed Endl et al. (2000). 
Briefly, DRV measurements were made with the iodine cell 
inserted in the light path which produces a superposition 

3 To calculate the RV secular acceleration we combined the 
Hipparcos proper motion and distance with the ARV value by 
Nidever et al. (2001) transformed to the same epoch. 
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Fig. 2. Time series of the differen- 
tial RV data of Barnard's star plotted 
vs. Barycentric Julian Day (bottom) 
and the Date (top). The solid line repre- 
sents the predicted RV secular acceler- 
ation of 4.499 m s _1 yr -1 . The dashed 
line is a linear fit to the data with a 
slope of 2.97 ± 0.51 m s _1 yr _1 . 



of a forest of dense iodine absorption lines onto the stellar 
spectrum. A higher-resolution (R = 700 000 - 1 000 000) 
template spectrum of the iodine cell obtained with the 
McMath Fourier Transform Spectrometer (FTS) on Kitt 
Peak was used to reconstruct the instrumental profile (IP) 
from observations of rapidly rotating B-stars through the 
iodine cell. This was required in one night for Maximum 
Entropy deconvolution (with the IP) of a pure stellar spec- 
trum (taken without the iodine cell) in order to obtain 
a stellar template spectrum with a resolution matching 
that of the FTS iodine spectrum. Finally both template 
spectra were combined to model the regular observations 
(star+iodinc), again employing IP reconstruction. From 
the Doppler shift of these model spectra with respect to 
the stellar template we obtain a time series of the DRV of 
the star which is largely corrected for instrumental insta- 
bilities. 

As the IP varies strongly within an Echelle order of 
UVES spectra as well as across the orders, the modelling 
was done in small spectral chunks of w 0.2 nm each of 
which yields an independent DRV measurement. DRV val- 
ues from 580 spectral segments (after rejection of out- 
liers) were then averaged to produce the final DRV value 
for the spectrum; the error of this mean value was taken 
as the measurement error. 

A total of 137 star+iodine spectra are available from 46 
nights between 21 April 2000 and 23 Sept. 2002. Triplets of 
exposures were obtained each night within about 16 min, 
except for one night where only two spectra could be se- 
cured. Individual exposure times were 250 — 285 s yield- 
ing a typical S/N per pixel of 66. On average our error of 
the individual RV measurement is 2.65 m s -1 demonstrat- 
ing the high precision achievable with UVES on this faint 
star (V = 9.54 mag). All RV data were corrected to the 
solar system barycenter using the JPL ephemeris DE200 
(Standish 1990) for the flux-weighted temporal midpoint 



of the exposure as given by the UVES exposuremeter. For 
each epoch of observation proper motion corrected stellar 
coordinates were used. 

A note on error estimation is in order at this point. 
As a multi-user instrument UVES in service mode un- 
dergoes frequent setup changes (central wavelength, slit 
width, image slicers, etc.). This can lead to systematic 
RV shifts of purely instrumental nature that may not be 
fully accounted for by our IP modelling approach. The IP 
will be more and more affected by setup changes when 
their number increases. Conversely, within a given night 
or on a time scale of a few nights higher instrumental sta- 
bility can be expected. Errors estimated from the scatter 
in a single night are most likely not representative for the 
long-term error. On the other hand, errors estimated over 
time scales of a few nights may already be affected by real 
variability of the star. 

For this reason we chose to perform all subsequent 
analysis of our data on nightly mean values of the DRV 
(46 data points) in order to avoid an overcstimation of 
the significances of possible signals. For the same reason 
we chose not to use the propagation of the errors of the 
individual nightly DRV values to obtain the error of the 
nightly mean DRV (which on average for all nights would 
be 1.54 m s _1 ); neither did we choose the RMS of the 
individual nightly DRV values around the mean for this 
purpose (which on average would give 1.39 m s" 1 ). Instead 
we chose the mean of the individual errors to represent the 
error of the nightly mean DRV (on average 2.65 m s _1 ). 
This choice may appear a bit conservative and arbitrary, 
but we estimate that it leads to more realistic values for 
the nightly errors as it is to a large extent the imperfec- 
tion of our IP reconstruction that is reflected in the error 
of the individual measurement and should therefore also 
be reflected in the nightly error. 
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Fig. 3. A selection of nightly mean 
spectra of Barnard's star near H a 6563 
and Cal 6573. Upper panel: three spec- 
tra are shown vertically offset. They 
correspond to the minimum H a index 
(lower spectrum), the maximum "qui- 
escent" H a index (middle spectrum off- 
set by 0.5 flux units), and the flare 
event (upper spectrum offset by 1.0 
flux units). Lower panels: Zooms into 
the upper panel on H a (left) and Cal 
(right). While a filling- in of the H a 
line by emission can be seen, the Cal 
line is constant. Horizontal bars indi- 
cate the spectral regions over which av- 
erage fluxes are determined for the line 
indices according to Eq. (5). 



Table 2. Models for the UVES data of Barnard's star: con- 
stant, predicted secular acceleration (S.A.), and linear fit. For 
each model we list the slope, the RMS scatter of the data 
around the model, the \ of the fit, the degrees of freedom 
(dof), and the associated probability. The RMS values should 
be compared with the mean internal measurement error of 
2.65 ms" 1 . The last line corresponds to a model based on the 
predicted secular acceleration together with a linear correlation 
with the H a line strength index defined in Sect. 5.2. 



Model 


Slope 


RMS 


x 2 


dof 


P 




[ms _1 yr _1 ] 


[ms" 1 ] 








Const. a 


(0.0) 


3.93 


97.2 


45 


1.3 ■ 10" 5 


S.A. a 


4.499 


3.38 


72.3 


45 


0.0065 


Lin. fit a 


2.97 ±0.51 


3.17 


63.4 


44 


0.029 


Lin. fit b 


5.15 ±0.89 


2.70 


29.0 


28 


0.41 


S.A.±H a c 


(4.499) 


2.94 


54.4 


43 


0.12 



a Data from all 46 nights. 

b Data from the first 2 years only (first 30 nights). 

c Data from all nights excluding flare event (45 nights). 



5. Results 

5.1. RV secular acceleration 

Fig. 2 shows the time series of our DRV data of Barnard's 
star. Each data point represents the nightly average of 
the DRV. The predicted secular acceleration of the RV 
(Sect. 3) is overplotted as a solid line (zero point matched). 
The dotted line represents a linear fit to the data. 

The best- fit linear trend yields a slope of 2.97 ± 
0.51 m s -1 yr -1 somewhat smaller (by 3. Oct) than the pre- 
dicted secular acceleration of 4.499 m s _1 yr -1 . Clearly, a 
constant model is rejected as it would be 5.8ct away from 
the fitted slope. For each model the slope and RMS scatter 



of the data around it are listed in Table 2 together with 
the x 2 , degrees of freedom, and associated probability. 

The relatively small formal probability of the secular 
acceleration model (0.65%) indicates that additional vari- 
ability is present in our data. The linear fit attempts to ac- 
count for part of this variability. If only the data of the first 
2 years are used for a linear fit (BJD < 2 452 000, i.e. the 
first 30 nights) the resulting slope is 5.15±0.98 m s -1 yr -1 
which agrees within 0.95ct with the predicted secular ac- 
celeration (not shown in Fig. 2). Around BJD = 2 452 500 
we had intensified the observations in an attempt to ob- 
tain a better sampling of a suspected minimum in the RV 
residuals from the secular acceleration model. The pres- 
ence of this minimum had been suggested by period analy- 
sis and the minimum was clearly found, but it added more 
weight to negative residual DRV values thereby reducing 
the slope of the linear fit. 

5.2. Correlation between H a line strength and RVs 

In order to find out whether stellar activity is the source 
of, or a contributor to, the observed RV variability we ex- 
amined the filling-in by chromospheric emission of the H Q 
absorption line in an attempt to find a correlation with 
the DRV data. H Q (6562.808 A) is the only activity indi- 
cator in the wavelength regime encompassed by our UVES 
spectra. Filling-in of the core of this absorption line is an 
indicator for chromospheric plage, i.e. active regions where 
magnetic fields are important for the convective proper- 
ties of the star and where star spots may also influence 
the shape of photosphcric absorption lines. 

M dwarf spectra are very rich in absorption lines which 
blend and do not show a continuum between them which 
makes it difficult to obtain precise absolute line fluxes. 
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Fig. 4. DRV (after subtraction of the secular acceleration) 
vs. spectral line index (see text). Circles with RV error bars: 
DRV vs. H a index. Crosses: DRV vs. Cal index. The line rep- 
resents a linear fit to the DRV vs. H a index relation. The right- 
most data point correpsonds to the flare event. 
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Fig. 5. Comparison of DRV time series and time series of H Q 
indices. Top: Time series of H a indices. Center: DRV time series 
(secular acceleration subtracted); RMS= 3.39 m s _1 (exclud- 
ing flare event). Bottom: Time series of DRV values corrected 
for the linear correlation with the H a index displayed in Fig. 4; 
RMS= 2.94 m s" 1 (excluding flare). 



Since, however, only relative fluxes are of interest here we 
define a suitable line index by 



0.5* (Fi +F 2 ) ' 

where F Q is the mean spectral flux in a selected RV in- 
terval around the line to be studied. For the H Q line we 
chose the RV interval [—15.5, +15.5 km s -1 ] centered on 
the core of the line. F± and F2 are the mean fluxes in 
two reference bandpasses on either side of H a (assumed 
to exhibit only negligible variability with the appear- 
ance of active regions) for which we chose the intervals 
[-700,-300 km s _1 ] and [+600, +1000 km s _1 ], respec- 
tively. These velocities are corrected for the barycentric 
Earth motion (cf. Sect. 3) and for the ARV of Barnard's 
star (Sect. 2). 

For comparison we evaluated the line index for a 
nearby Cal line at 6572.795 A with low excitation po- 
tential that is very strong in M dwarfs and whose line 
strength is not expected to have a significant dependance 
on stellar activity. For this line we chose F Q over the 
RV range [+441.5, +472.5 km s _1 ] (again relative to H a ). 
Should our data reduction have produced systematic flux 
offsets, e.g. due to imperfect subtraction of the Echellc 
background, we should see them in the Cal line. 

Fig. 3 shows three of our (nightly mean) spectra of 
Barnard's star near the H a 6563 and Cal 6573 lines. They 
correspond to the spectrum with the smallest H Q index 
(labelled "Min."), to the spectrum with the largest "qui- 
escent" H Q index (labelled "Max.") and to the spectrum 
of a flare event with an H Q index that is almost 12cr larger 
than the mean index of the remainder of the data (see 
also Fig. 4). Horizontal bars in Fig. 3 indicate the spec- 
tral regions over which the mean fluxes for Eq. (5) were 
determined. 



The variation of the H a line strength can be seen as 
filling- in by emission, whereas the Cal line is constant. 
Interestingly, the smallest H Q index was observed in the 
night following the flare event. It appears that either the 
flare occured on a quite inactive (visible) hemisphere or 
that the flare event has disrupted a major active region 
(just as solar flares cause a reconfiguration of the local 
magnetic fields). Since we have no data from the night be- 
fore the flare these possibilities must remain speculation. 

Fig. 4 shows a correlation of our DRV measurements 
(secular acceleration subtracted) with the H Q and Cal in- 
dices. Just as the DRV values also the line indices were 
determined in the individual spectra of each night and 
then binned into nightly averages. Again, the Cal line is 
seen constant and independent from the DRV, but the H Q 
line strength is variable and anti-correlated with the DRV. 
The flare event is seen only in Hq,, but not in the Cal in- 
dex. The relative scatter of the line indices is 0.54% for Cal 
and 7.6% for H a (or 3.8% when the flare is excluded). The 
correlation coefficient is —0.498 (flare event excluded). A 
linear fit to the relation of DRV vs. H a index, I(H a ), is 
given by DRV = (43+10) m S - 1 -(60±15) m s" 1 7(H Q ) . 
It reduces the total RMS of the DRV data (flare excluded) 
from 3.39 m s _1 to 2.94 m s _1 . The probability of x 2 for 
this linear fit is p = 0.12 rendering this model an ac- 
ceptable description of the DRV data (see Table 2 for all 
numerical values of the fit). 

A technical note: near H Q our spectra still con- 
tain weak iodine absorption lines from the iodine cell 
(cf. Sect. 4). The deepest iodine line in the wavelength 
regime relevant for our line indices has a depth of 9.3% 
relative to the continuum; about half of this value is more 
typical for the average line. Due to the barycentric mo- 
tion of the Earth, for which the chosen RV intervals are 
corrected (with correction extrema of —26.4 km s _1 and 
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Fig. 6. Lomb-Scargle periodograms for 
our data of Barnard's star. Top panel: 
H a index. Second panel from top: RV 
data (secular acceleration subtracted). 
Third panel: RV data additionally cor- 
rected for the correlation with the H a 
index. Bottom panel: window function. 
In each case the flare event was ex- 
cluded from the data. Horizontal dotted 
lines indicate levels of the FAP (1% and 
5% in the upper two panels and 10% 
and 50% in the third). The vertical dot- 
ted line markes the supposed 130.4 d ro- 
tation period by Benedict et al. (1998). 



+26.2 km s -1 ) variable parts of the iodine spectrum con- 
tribute to the flux values in Eq. (5). Calculating the H Q 
and Cal indices for a pure iodine spectrum (a Hatfield ex- 
posure taken through the iodine cell) we obtain an RMS 
variation of 0.20% and 0.17% for the H Q and Cal index, re- 
spectively, over the range of applied barycentric correction 
velocities. This intrinsic uncertainty is much smaller than 
the observed scatter of the line indices for our spectra of 
Barnard's star (3.8% and 0.54%, respectively). Therefore, 
our neglection of the contribution of the weak iodine lines 
to these spectra does not introduce a major error in the 
determination of the line indices. 

In Fig. 5 we show a comparison of the time series of our 
DRV data (secular acceleration subtracted) with the time 
series of the H Q index. Also included in Fig. 5 is the time 
scries of the DRV data corrected for the variation with 
H Q , i.e. after subtracting the linear fit shown in Fig. 4. 

Even if the total scatter of the DRV data is reduced 
when the correlation with the H Q index is applied for a 
correction, the basic variability pattern appears largely 
unchanged. Formally, the probability of x 2 provides an 
acceptable model, but we cannot exclude the possibility 
that our error estimates were too conservative (cf. Sect. 4). 
We should also note that we do not expect a simple linear 
correlation to be sufficient for a complete correction of the 
RV data for the effects of stellar activity. 

5.3. Period analysis 

We searched for periodic signals in the H Q index, in the 
RV residuals from the secular acceleration model, and in 
the RV data after additional correction for the linear cor- 
relation with the H Q index. We used the periodogram by 
Lomb (1976) and Scargle (1982) as well as a routine that 
employs x 2 fitting of sine waves. Contrary to the Lomb- 
Scargle method the sine fitting routine can adapt the RV 



Table 3. Period analysis results for the different data sets: 
H a index, RV after subtraction of the secular acceleration, RV 
additionally corrected for the linear correlation with the H a 
index, and window function. For the data sets that exclude the 
flare event the period values P of the four highest periodogram 
peaks as well as the pertaining FAP are listed. Only for the RV 
data did we also evaluate period and FAP values with the flare 
event included. Similarity of a period value with an integer 
fraction of the photometric period by Benedict et al. (1998) is 
indicated in the last column. 



Data set 


P[d] 


FAP 


P [d] FAP 


Comment 




flare excluded 


including flare 




H a index 


44.5 


0.82% 




« ^ph/3 




30.4 


1.30% 








28.2 


1.54% 








50.4 


2.06% 






RV 


32.7 


0.56% 


32.7 1.56% 


w Pph/4 




45.3 


1.05% 


45.3 4.27% 


» Pph/3 








19.1 5.06% 






36.1 


4.01% 


36.0 8.26% 






51.1 


5.49% 






RV— H a 


32.5 


36.81% 




« Pph/4 




40.0 


45.21% 








19.1 


75.63% 








440 


77.96% 






Window 


380 






w 1 yr 




190 






« 1/2 yr 




52.5 










46.2 









zero point offset and also take into account data errors 
(cf. Cumming ct al. 1999). The results of both methods 
do not differ substantially so that we only present the re- 
sults from the Lomb-Scargle periodogram. Fig. 6 shows 
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Fig. 7. Upper limits to the RV semi-amplitude K (for circular 
orbits) as a function of orbital period P. For each period value 
these limits correspond to the semi-amplitude at which the 
artificial signal was found at the FAP < 1% level for all 12 test 
phases. Two vertical dotted lines at periods of 5.75 and 21.5 d 
indicate the range of the habitable zone around Barnard's star 
according to Kastings et al. (1991). Periods of 0.5 and 1 yr are 
also marked by vertical dotted lines; due to the structure of 
the data sampling relatively high limiting amplitudes result at 
these periods. 



periodograms of the H Q index, the RV residuals, and the 
RV data corrected for the correlation with H Q ; the window 
function is also shown. 

We searched the period range from 2 d (our minimum 
temporal resolution) to 2.42 yr (our total time baseline). 
False alarm probabilities (FAP) for the periodogram peaks 
were determined with a bootstrap randomization scheme 
(e.g. Kiirster et al. 1997). 1% and 5% levels of the FAP 
are indicated in Fig. 6 for the periodograms of the H Q and 
RV residuals. After correction with the H Q correlation the 
power of the RV data was found considerably reduced; we 
display 10% and 50% FAP levels for this data set. 

Table 3 lists the periods and FAPs of the four highest 
periodogram peaks for each of the three different data sets 
excluding the flare event as well as for the RV data includ- 
ing the flare event (pertaining periodogram not shown in 
Fig. 6). In none of the periodograms is there much power 
near the supposed photometric rotation period, 130.4 d, 
by Benedict et al. (1998) which is also indicated in Fig. 6 
as a vertical dotted line. Nevertheless, the highest peaks 
tend to coincide with 1/3 or 1/4 of the photometric period 
value which may hint at a common origin. 

5.4. Upper limits to companion masses 

Faced with the result that the observed RV variability 
appears to be attributable to stellar activity to a large 
extent and yields no evidence for orbiting companions we 
studied which types of companions can be excluded, i.e. we 
determined upper limits to the projected mass msmi of 
orbiting companions (with i the orbital inclination). To 



Fig. 8. Upper limits to the mass m of planetary companions 
in circular orbits around Barnard's star plotted vs. star-planet 
separation a (bottom) and orbital period P (top). Long-dashed 
line labelled "msini": 99% confidence limits to the projected 
companion mass. Solid line labelled "m": limits to the true 
companion mass based additionally on the 99% confidence 
range for the inclination i; i.e. their combined confidence is 
98%. Short-dashed line labelled "AL": astrometric limits to m 
from Benedict et al. (1999). Horizontal thin dotted lines indi- 
cate the masses of Jupiter, Saturn, Neptune, and Uranus. As 
in Fig. 7 the habitable zone according to Kasting et al. (1991) 
lies between the vertical thin dotted lines. 



this end we determined the amplitudes of RV signals that 
we would have detected with 99% confidence. 

Again bootstrap randomization was employed. In the 
analysis artificial RV signals (see next paragraph for de- 
tails) are added to the observed RV data which act as a 
noise term. Since on the one hand secular acceleration is 
a well understood effect, but on the other hand the lin- 
ear correlation with the H Q is probably not a complete 
description of the effect that activity exerts on RV mea- 
surements, we decided to use as the noise term the secular 
acceleration subtracted RV data including the flare event 
and uncorrected for the H Q correlation. This data set has 
a higher rms scatter than the corrected one (see Table 2), 
i.e. it adds more noise to the test signals and yields more 
conservative upper limits. Note also that the highest peak 
in the periodogram of this data set itself (no signal added) 
appeared at FAP = 1.56% (fifth column in Table 3) and 
would therefore not have passed the criterion of 99% con- 
fidence that we adopt here. 

Sinusoids were used for the artificial signals render- 
ing our analysis strictly valid only for circular orbits. We 
used 29 different period values encompassing the accessi- 
ble period range 2d — 2.42 yr. At each period 12 equally 
spaced phases were probed, and the amplitude of the sine 
wave was varied. The amplitude at which for a given in- 
put period the signal showed up at all 12 phases at the 
FAP < 1% level was taken as the minimum detectable 
amplitude. We did not require that the signal appear at 
the same period in the power spectrum (due to spectral 
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leakage it not always did), but that it cause power with a 
chance probability < 1% somewhere in the search range. 

Fig. 7 shows the resulting limiting RV semi-amplitudes 
as a function of orbital period. The median (mean) de- 
tectable RV semi-amplitude in the whole period search 
range is 5.88 (6.20) m s _1 ; in the habitable zone it is 
4.73 (5.07) His" 1 . These limits are smaller than any of 
the amplitudes found so far for extrasolar planets. 

The amplitude limits were then converted to upper 
limits on the projected mass mshn of planetary compan- 
ions in circular orbits. A stellar mass of 0.16 M Q was 
adopted from Henry et al. (1999). The limits on msinz 
are shown as a function of period and star-planet sepa- 
ration in Fig. 8 (long-dashed line). To obtain limits to 
the full companion mass m we used the probability (for 
random orientation of orbits) that the orbital inclination 
i exceeds some angle 9, given by p(i > 9) = cos(#). In 
99% of the cases the true mass does not exceed the pro- 
jected mass by more than a factor 7.088. The solid line 
in Fig. 8 shows these mass limits which have a combined 
probability (that of signal detection and that of the in- 
clination range ) of 0.99 2 w 98%. The short-dashed line 
labelled "AL" shows astrometric mass limits by Benedict 
et al.(1999) which are 99% confidence limits, but allow for 
a 5% miss rate for astrometric signals. 

As seen in Fig. 8 the limits on m sin i increase roughly 
from a few Earth masses at the smallest separations to 
i=a 30 M Earth at larger separations. At periods of 1/2 yr and 
1 yr we obtain somewhat higher limits because of lower 
detection capabilities due to the data sampling window. 

6. Discussion 

6.1. RV secular acceleration 

We clearly find a positive trend in our DRV data of 
Barnard's star and interpret it as mainly due to the secular 
acceleration of the RV. The formal discrepancy between a 
linear fit and the predicted secular acceleration value can 
be attributed to stellar acticity, in particular to incom- 
plete sampling of the time scale(s) on which the activity 
varies. 

To our knowledge this is the first reported observation 
of stellar RV secular acceleration. With time we expect 
the observed trend to conicidc more with the predicted 
acceleration as activity effects cancel out more. 

6.2. The activity-related RV signal: Evidence for 
convective redshift? 

Active regions, which in our data reveal their appearance 
through the filling-in of the H a line core by emission, can 
affect the apparent stellar RV as measured in lines of pho- 
tospheric origin by the effects of star spots or chromo- 
sphcric plage regions. Compared to the inactive photo- 
sphere both spots and plage produce a local change of 
brightness and spectrum and also inhibit stellar convec- 
tion locally by magnetic fields. 



The primary effect of dark star spots at different lo- 
cations on the visible stellar disk is a reduction of the 
observable flux at the instantaneous RV of the spot(s), as 
spot temperatures T s are considerably smaller than pho- 
tospheric effective temperatures T e g and (T s /T e g) 4 <C 1. 
Therefore, the effects that spots have a different spec- 
trum and that convection is suppressed inside the spot 
have small contributions to the integrated stellar spec- 
trum compared to the general flux reduction in a small 
RV regime. If there is sufficient rotational broadening of 
the absorption lines, then this will lead to a change in the 
center-of-gravity of the line thus mimicking an RV change. 
The effect is more pronounced when spots are displaced 
from the limb of the visible stellar disk where their in- 
stantaneous RV (relative to the RV of the star as a whole) 
is high. The effect disappears when spots are on the sub- 
observer meridian of the star where the relative RV is zero 
(or when they are on the invisible hemisphere of the star). 

In a model that assumes that photospheric spots are 
surrounded by chromospheric plage the following picture 
emerges. Minimum RV shift, i.e. most active regions hid- 
den from view or most active regions visible and con- 
centrated towards the sub-observer meridian, should be 
accompanied by either small or large amounts of the 
filling-on of the H Q line core by emission. Maximum RV 
shift, i.e. active regions concentrated on either side of the 
sub-observer meridian, should be seen when intermediate 
to large amounts of H Q filling-in are observed (see also 
Paulson et al. 2002; Saar & Fischer 2000). 

Consider the simple case of a single small spot on an 
otherwise featureless rotating star. Let A s be the spot area 
and let <f) be the rotation angle of the spot from the sub- 
observer meridian. The RV perturbation Av r then scales 
as Avr((j)) oc —A s f(<f))g(<f))v sini, where ^4 S is the spot area, 
f{4>) is a function describing the visibility and foreshort- 
ening due to projection of the spot, hence A s f(<f>) is the 
projected spot area, g(4>) is a function describing the de- 
flection of the stellar RV caused by the spot as a function 
of stellar rotation phase 0, and wsini is the projected ro- 
tational velocity of the star. The minus sign indicates that 
the RV shift of the line as a whole is opposite to the in- 
stantaneous RV of the dark spot due to the flux reduction 
it produces. 

Both f{4>) and g(<j>) are periodic functions of <f) which 
possess special symmetries around = 0, i.e. the phase 
where the spot is on the sub-observer meridian. The pro- 
jection function f{<f>) is mirror symmetric with respect to 
this phase (similar to a cosine function) while the RV 
deflection g(<p) has a rotational symmetry around this 
phase (just like a sine function does). In addition g(<fi) 
has zero mean. 4 Due to these properties of f(<f>) and 
g{4>) the RV perturbation averaged over phase <j) vanishes, 
i.e. (Av T ) = 0. On time scales longer than the stellar ro- 
tation period the mean spot related RV perturbation is 

4 For the RV deflection to be a fully rotationally symmetric 
function with zero mean the absorption line(s) under study 
must have symmetric profiles which is a good approximation. 
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therefore uncorrelated with the visibility of spots and spot 
coverage as indirectly evidenced by H Q emission from the 
surrounding plage. It is only the scatter in the instanta- 
neous RV deflection which correlates with spot visibility 
and area, but not the mean value. 

In a plot of DRV vs. H Q index one should thus ex- 
pect for a single active region (small spot plus surrounding 
plage) a closed curve with a shape that resembles a "re- 
clining raindrop" ; 5 if several active regions are present 
the plot contains a superposition of curves of this shape. 
If plage regions exist without underlying spots, the H Q in- 
dex will be further enhanced when they are visible without 
any spot-related effect on the RV, thus blurring the plot 
along its horizontal axis. In any case, if spots are responsi- 
ble for the RV variability the correlation coefficient of the 
RV vs. H a relation should be close to zero. This is not ev- 
ident in Fig. 4 so that star spots cannot be the dominant 
contributor to the observed correlation. 

Estimates of the maximum spot contribution to Av r 
are consistent with this. Saar & Donahue (1998) show 
that the spot-induced RV amplitude (in m s _1 ) is K s w 
6.5/ s °- 9 usini, where f s is the differential spot filling fac- 
tor (in %) and vsini is in km s _1 (see also Hatzes 2002). 
Since the stellar radius is roughly R ?a 0.2 Rq and the 
rotation period is P ro t > 32.5 d (assuming that the pe- 
riod of highest power in the periodogram of the RV data 
in Table 2 is either the rotation period or a fraction 
thereof), one obtains vsini < 0.31 km s^ 1 . If the spot- 
related brightness variation of the star is still similar to 
that found by Benedict et al. (1998), i.e. AV w 0.01, im- 
plying a change in spot filling factor of / s w 0.92%, we get 
K s < 1.9 m s _1 , and the associated spot-induced scatter 
<7 S w K/ \pi < 1.3 m s _1 , far lower than the variation seen 
here. 6 This is further indication, under the assumptions 
we have made, that spots are not the main contributor for 
the RV variability in Barnard's star. 

The effects of bright chromospheric plage regions on 
the integrated spectrum of the star are quite different from 
those of spots. In the visual the plage-photosphere bright- 
ness contrast is small and does not exert a strong effect 
on photospheric absorption line shapes. Also such an ef- 
fect would cause similar cyclic RV shifts as spots (with an 
inverted sign) and would not lead to a correlation with the 
RV. However, the brightness of plage lets them contribute 
their own spectrum to the overall stellar spectrum. For 
the same reason the suppression of convective flows in the 
plage magnetic fields can lead to observable effects. 

In solar-type stars one usually observes the plage spec- 
trum in the form of enhanced emission in those lines in 
which chromospheric emission is most readily seen, such 

5 The shape of a "reclining raindrop" is obtained when the 
popular picture (not the physical one) of a raindrop is rotated 
counterclockwise by 90°. 

6 Subtracting in quadrature the total scatter of the RV data 
after subtraction of the secular acceleration, 3.38 m s _1 , and 
the average internal error of 2.65 m s _1 (conservative estimate; 
cf. Sect. 4) we obtain an excess scatter of 2.10 m s _1 . 



as Call H+K or the Balmer lines. This emission is usually 
blueshifted with respect to the photospheric line compo- 
nent and its strength increases with activity level thereby 
changing the degree of line asymmetry. If a similar ef- 
fect also ocurrs in other photospheric lines (perhaps to a 
smaller degree), then it could affect the measured RV. So 
far this has not yet been conclusively demonstrated. 

However, in our H a line profiles of Barnard's star we 
find no evidence for a distinct behaviour of the blue or 
red parts of the line. To this end we measured the ratio of 
the mean fluxes in the RV intervals [—35,0 km s _1 ] and 
[0, +35 km s _1 ] around H Q (cf. Fig. 3) and found it con- 
stant within 1.9% rms and uncorrelated with both our RV 
data and the H Q index. If a variable differential emission 
enhancement between the blue and red parts of the H Q 
line of Barnard's star exists, then its degree is below our 
measurement precision. 

Clearly, we cannot draw direct conclusions from this 
finding on the contribution of the chromospheric plage 
spectrum to the absorption lines relevant for our RV mea- 
surements. Nevertheless a tentative assumption that such 
a contribution is small does not seem unreasonable. 

Then what about the suppression of the local convec- 
tive velocity pattern by magnetic plage? Local suppression 
of convection alters the center-of-gravity of photospheric 
absorption lines (see Cavallini et al. 1985 for the solar 
case). If the fundamental pattern is convective blueshift, 
as seen in main sequence stars of types G or K, then any 
local suppression of the pattern produces a net redshift of 
the lines when active regions appear. 7 

To repeat our basic assumptions, we consider the 
case in which that part of the spectrum that is rele- 
vant for our RV measurements is similar when emitted 
either from plage or from the inactive photosphere, both 
in continuum brightness and absorption line strengths and 
shapes. Then we find that for each individual plage region 
the plage induced RV perturbation scales as Av T (<f)) oc 
—A p v cs ((f))f(<p), where A p is the surface area of the plage, 
v cs (<p) is the convective blueshift in the absence of plage, 
and f{<j>) describes the visibility and foreshortening of the 
plage area. Averaging over the rotation phase </> we obtain 
(Ai> r ) cx -A p (v cs (4>)f(4>)}. Since neither v cs (4>) nor f(4>) 
have zero mean the plage induced RV redshift does not 
vanish and increases with plage area. If more than one 
plage region is visible, their contributions are additive. 

In this case there should be a clearer correlation be- 
tween the RV and the total plage coverage as observable 
in the H Q index. This seems to be similar to what we see 
in Fig. 4, however, with opposite sign. Our figure reveals 

7 For reasons of the conservation of the energy flux sup- 
pressed convection in magnetically active regions should in 
principle be compensated by an increase outside of them. To 
our knowledge searches for "bright rings" around sunspots, 
where the "missing energy" was expected to emerge, have not 
been successful. In analogy with the Sun we implicitely assume 
that the reduced convective line shift due to the local suppres- 
sion of convection in plage dominates over the effect of possibly 
increased convection outside. 
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a net blueshift when the coverage of the star with plage 
regions increases which can be explained, if the convection 
properties in Barnard's star produce convective redshift as 
the fundamental pattern, i.e. in regions not dominated by 
magnetic fields. If this fundamental pattern is suppressed, 
a net blueshift will be the result, as observed. 

The possibility that M dwarfs have convective redshifts 
may extend a trend (by way of extrapolation) seen in re- 
cent work. Gray (1982) found that bisectors in dwarfs be- 
come less curved ( "C shaped" ) and more vertical (see also 
Saar & Bruning 1990), and also less blueshifted with de- 
creasing effective temperature T c ff . Although Gray's veloc- 
ity zero-point was somewhat uncertain, his data possibly 
hints that by mid-K main sequence stars, the net convec- 
tive blueshift may actually change sign to a convective 
redshift (see also Dravins 1999). 

If convective redshift prevails in Barnard's star, then 
the majority of the photosphcric absorption lines relevant 
for our RV measurements seems to form in regions of con- 
vective overshoot. Another possibility is that the details 
of the convection pattern, i.e. the interplay of the involved 
flow velocities, contrasts and surface areas of granules 
(where material rises) and inter-granular lanes (where the 
gas flows down) cause a net redshift of the relevant lines. 

Theory also supports the idea of a decreasing convec- 
tive blueshift with decreasing T c g (Dravins & Nordlund 
1990, Asplund et al. 2000). Although analogous line pro- 
file calculations for 3-D hydrodynamic M dwarf mod- 
els have not yet been published, the models themselves 
(Ludwig et al. 2002) show low contrast between upflow- 
ing and downflowing regions, C(dM) w 1.1%, and low 
RMS vertical velocities of up to (7j_(dM) w 0.24 km s _1 . 
Analogous solar models (see also Ludwig et al. 2002) 
yield contrasts C(O) «16% and RMS vertical velocities 
<7j_(0) w 2.6 km s _1 . Adopting these values permits us to 
make a rough estimate of the expected RV shift due to 
convective motions, v CSl if we assume that it scales with 
<7j_. Then we can expect (including the sign change) 



Vcs(dM) 



<r±(dM) (1 - C(dM)) 

'~^W (i-c(©)) 



(6) 



With w cs (0) w -300 m s" 1 (e.g. Dravins 1999), this im- 
plies a convective redshift of w cs (dM) w 33 m s _1 in the 
absence of plage. If plage inhibit the convection at vari- 
able amounts as they appear, grow, shrink, disappear, or 
rotate out of view, the possible variation of the convective 
shift cannot exceed Av cs (dM) w 33 m s" 1 . 

This value is consistent with our observed (smaller) 



range of RV variation of 6.9 ± 1.7 m s 



obtained from 



the linear fit to the relation between DRV and H Q index 
(Fig. 4). We thus estimate that the variation in the surface 
filling factor of visible plage is approximately (21 ± 5)%. 

6.3. Planetary companions 

Our detection capabilities (Figs. 7 and 8) with UVES are 
currently leading among the ongoing precision RV surveys. 



To our knowledge the mass upper limits determined in 
Sect. 5.4 for planetary companions in circular orbits are 
the lowest published so far. 

For circular orbits we find no planets with to sin z > 

0. 12 Mj U pitor and m > 0.86 Mj up i tcr over the range of 
separations accessible to our RV data (0.017 — 0.98 AU) 
corresponding to periods of 2—885 d. In the habitable zone 
around Barnard's star (separations of 0.034 — 0.082 AU 
and periods of 5.75 — 21.5 d), we can even exclude planets 
with to sin i > 7.5 MEarth and to > 3.1 MNcptunc- 

We cannot fully exclude the possibility that we 
may have missed planets with higher masses in high- 
eccentricity orbits. Note, however, the astrometric limit 
in Fig. 8 (from Benedict et al. 1999). Taking into account 
eccentric orbits in our simulations would make the calcu- 
lations prohibitively expensive, because of two additional 
free parameters (eccentricity and longitude of periastron). 
For further discussion see Endl et al. (2002). 

7. Conclusions 

1. We have clearly discovered a linear trend in the RVs 
of Barnard's star which we attribute to the RV secular 
acceleration of the star. 

2. The residuals from the secular acceleration are corre- 
lated with stellar activity to a large extent. 

3. From the sign of this correlation we propose that the 
fundamental (magnetic field free) convection pattern 
in this M dwarf is convective redshift as opposed to 
the blueshift seen in earlier spectral type stars. 

4. A rough calculation based on the convective redshift 
assumption leads us to estimate that the variation of 
the visible plage coverage is about 20%. 

5. From our RV data we infer the so far lowest pub- 
lished mass upper limits for planets in cicrular or- 
bits around their host star. We exclude planets with 
to sin i > 0.12 Mj up i te r and to > 0.86 Mj up i tor over 
separations of 0.017 — 0.98 AU. In the habitable zone 
around Barnard's star, i.e. 0.034 - 0.082 AU, our RV 
data exclude planets with to sin? > 7.5 MEarth and 
to > 3.1 M Nopt line- 
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